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Abstract: The flux of neutrinos and photons originating from cosmic ray inter- 
actions with the interstellar medium in the galaxy is calculated based on current 
models for high energy particle interactions. The contribution from a possible dark 
matter halo of the galaxy is considered. The photon flux gets a non-trivial attenua- 
tion due to interactions with the cosmic background radiation. The neutrino fluxes 
are compared with those originating from the Earth's atmosphere as well as from 

s active galactic nuclei. 

^ ' 

o 

! 1 Introduction 

Cosmic rays of galactic and extra-galactic origin will interact in high energy collisions with the 
interstellar medium of our galaxy and produce secondary particles (l|, ||, ||. These are mainly 
mesons that decay and give rise to a flux of muons, neutrinos and photons. The very low 
,^Zh • density of the interstellar medium imply that the interaction lengths of the secondary particles 

is long compared to their decay length, such that the mesons will decay before loosing energy 
in secondary interactions. This is also the case for the muons which will decay giving neutrinos. 
This is the opposite of the situation for cosmic ray particles interacting in the Earth's atmosphere 
Q] , where meson typically loose energy in interactions before decaying. The fluxes of high energy 
neutrinos and photons from the interstellar medium could therefore be larger than those from 
the atmosphere, although the initial production rate of mesons is smaller. 

A measurement of these fluxes could potentially give valuable information about the distribu- 
tion of matter and cosmic rays in the galaxy, which could be of great importance in determining 
the origin of the cosmic rays. Understanding the flux from the disc of the Milky Way could also 
be the starting point in a search for baryonic dark matter in a spherical halo around the Milky 
Way§|. 

In addition, these fluxes from the interstellar medium constitute a background in searches 
for other, more spectacular cosmic sources and must therefore be known in order to extract the 
desired signal. For example, there is much interest in neutrinos from Active Galactic Nuclei 
(AGN) (see for example || and references therein). 

In this paper we present a realistic calculation of the neutrino and photon fluxes to be 
expected from cosmic ray particle interactions in the interstellar medium. These fluxes are 
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derived from complete events obtained by detailed Monte Carlo simulations based on state-of- 
the-art models for high energy particle collisions. In section 2 we present the model including 
the cosmic ray energy spectrum, the interstellar matter distribution in the Milky Way and the 
model for particle production in high energy elementary particle collisions. The resulting fluxes 
are presented in section 3, where also the attenuation of the photon flux due to interactions with 
the cosmic background radiation is demonstrated to give a significant effect. Section 4 considers 
the mentioned dark matter halo and the fluxes that would arise from it. We end, in section 5, 
with a discussion and comparison with estimated neutrino fluxes from active galactic nuclei. 

2 The calculational method 

The fluxes of neutrinos and gammas that could be observed at the Earth is basically given by 
the production rate in the galaxy and the probability that the produced particle reaches the 
Earth. The fluxes can be expressed as an integral along the line of sight from the Earth (r = 0) 
to the edge of the galaxy (r = R) 

where the primary cosmic ray flux <p p (E',r) is folded with the interstellar matter density p(r) 
and the differential cross section for producing neutrinos or photons. An attenuation factor 
A(E, r) is also included to account for the loss of flux due to interactions with the interstellar 
medium and P(r) is the probability that a particle is directed toward the Earth. A(E, r) is 
calculated separately and folded with the non-attenuated fluxes. Its effect is important for the 
gamma flux were interactions with the cosmic background radiation producing electron-positron 
pairs occur. 

The integration over r can be greatly simplified since one may calculate the flux from a 'unit' 
column density and then scale with the appropriate integrated density in any given direction. 
This procedure is possible since the distance between a typical production point r and the Earth 
is large compared to the decay length for all particles. The energy integral can be obtained 
based on Monte Carlo simulations of high energy particle collisions giving complete final states 
including particle decays. For this purpose we have employed the Lund Monte Carlo programs 

1- 

2.1 The cosmic rays 

The flux of cosmic rays in the Milky Way is not well known. The flux measured at the Earth 
is found to be isotropic to a high degree, the anisotropy being ^5% 0. This isotropy can arise 
in two ways. First, the bulk of the flux is of extra-galactic origin with a small component from 
processes in the Milky Way, i.e. the flux is universal. Secondly, it is of galactic origin with only 
a component at the highest energies of extra-galactic origin. The flux at lower energy is in this 
case captured by the magnetic field of the galaxy and thereby appearing to be isotropic at the 
Earth. A flux of extra-galactic origin will not be significantly attenuated when passing through 
the galaxy, since only a few per mill of the cosmic ray particles will interact. This attenuation 
effect is smaller than other uncertainties concerning the primary flux, and the flux can be treated 
as homogeneous. 

Assuming that the Earth is not at a unique site in the galaxy, it follows that the flux is 
the same everywhere in the galaxy. Therefore, we assume the flux observed at the Earth to be 
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uniformly and isotropically distributed throughout the galaxy. The energy dependence of the 
flux can then be parameterized as M, E| 
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nucleons 
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1.7 E~ 2 - 7 £<5-10 6 GeV 
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The normalisation is here derived [10] from the directly measured primary spectrum using 
balloon-borne emulsion stacks in JACEE [11|. It agrees (within some 10%) with more indi- 
rectly derived spectra based on measured atmospheric muon fluxes |12]] , and is also compatible 
with the data discussed in ref. [13|. The cosmic ray composition is dominated by protons with 
only a smaller component of nuclei [10, fill]. 



2.2 The model of the Milky Way 

Given the uncertainties concerning the matter distribution in the galaxy, one can only make a 
fairly simple model for it. However, this should also be adequate for our purpose to estimate 
the fluxes. Such a model has been used by Domokos et al. but it seems to overestimate the 
column density for high galactic latitudes. Therefore, a model similar to the one in j|] will be 
used outside the plane of the Milky Way, where the density decreases exponentially with the 
distance to the galactic plane. 

We take the galaxy to be rotationally symmetric in the plane with a radius of 12 kpc and 
with a constant constant density of lnucleon/cm 3 in the plane. Outside the plane, the density 
decreases as 

p{h)=p e- h l h \ (3) 

where h is the height above (or below) the disc and ho = 0.26 kpc is the scale height at the galacto- 
centric distance of the Earth. The Earth is in the galactic plane at a distance of Roc = 8.5 kpc 
from the centre. 

Stars are not taken into account in this study, although they affect the fluxes in two ways. 
First, the cosmic ray particles can interact in the stellar material producing mesons. These will, 
however, often have secondary interactions due to the relatively higher density environment 
in the stars. They will thus typically be absorbed or degraded in energy before they decay 
and therefore relatively less important [pj . Furthermore, the neutrinos and gammas from these 
meson decays must traverse the star such that photons will be absorbed and also the neutrino 



flux be significantly attenuated. These aspects are treated in our study |14] of the neutrino 
fluxes arising from cosmic ray interactions in the Sun. Secondly, neutrino and gamma fluxes 
from the interstellar medium will be attenuated because of absorption in stars. However, this 
is a very small effect since the solid angle covered by stars is only AO/f^lO -10 (derived from 



2.3 The model for particle production 

To specify the energy spectra of secondaries in cosmic ray collisions with the nuclei of the 
interstellar medium a model for particle production is needed. With the cosmic rays being 
predominantly protons and the interstellar matter mainly hydrogen, and only some smaller 
fraction of heavier nuclei in both cases, the interactions producing the secondary particle fluxes 
are dominantly proton-proton collisions. Contributions of nuclear collisions can also be treated as 
nucleon-nucleon interactions, since the nuclear binding energies are negligible and other nuclear 
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effects have only little influence on the high energy secondary particles that are of interest to 
us. Therefore, the primary interactions can be taken as proton-proton collisions. 

The dominating source of neutrinos and gammas are then the decays of the light mesons (ir 
and K) and muons plus a small contribution from other heavier hadrons. Given the low density of 
the interstellar medium, secondary interactions are rare and will essentially not happen before 
unstable particles decay. Therefore, with all light hadrons decaying, the contribution from 
charmed and heavier particles will not be important as opposed to the case of interactions 
in the Earth's atmosphere @. The dominance of light hadrons (pions) will only disappear 
above their critical energy ^<^" ltlcal 10 21 GeV ]l]]), where their decay becomes less probable 
compared to their interaction. Heavy quark production is therefore not explicitly considered 
in the simulation, but their contribution is estimated based on the spectrum-weighted moment 
method 0] and shown to be negligible. 

The production of light hadrons (not containing heavy quarks) is dominantly through min- 
imum bias hadron-hadron collisions. The strong interaction mechanism is here of a soft non- 
perturbative nature that cannot be theoretically calculated from first principles, but must be 
modelled. In the successful Lund model |l5|] hadron production arise through the fragmenta- 
tion of colour string fields between partons scattered in semi-soft QCD interactions Q. The 
essentially one-dimensional colour field arising between separated colour charges is described by 
a one-dimensional flux tube whose dynamics is taken as that of a massless relativistic string. 
Quark-antiquark pairs are produced from the energy in the field through a quantum mechanical 
tunneling process. The string is thereby broken into smaller pieces with these new colour charges 
as end-points and, as the process is iterated, hadrons are formed. These obtain limited momenta 
transverse to the string (given by a Gaussian of a few hundred MeV width) but their longitu- 
dinal momentum may be large since it is given by a probability function in the fraction of the 
available energy-momentum in the string system taken by the hadron. All mesons and baryons 
in the basic multiplets may be produced and the subsequent decays are fully included. The 
iterative and stochastic nature of the process is the basis for the implementation of the model 
in the Jetset program [||. The parameters in the program were taken at their default values, 
except that also long-lived particles had to be treated as unstable and their decay simulated. 

A non- negligible contribution to the inclusive cross section is given by diffr active interactions. 
These are also modeled in Pythia || using cross sections from a well functioning Regge-based 
approach and simulating the diffractively produced final state using an adaptation of the Lund 
string model. These diffractive events are included in our simulations and contribute rather less 
than 10% to the final results. 

3 Results 

As mentioned, the production rate is only a function of the density since the cosmic ray flux 
is assumed to be position independent within the galaxy. The energy integral of Eq. (|l|) has 
been calculated for a 'unit' column density of 1 GeV/cm 3 over a distance of 1 kpc (~ 5 mg/cm 2 ). 
This was performed by choosing energies between 10 2 and 10 10 GeV from a flat distribution 
in logioE, and assigning a weight to account for the detailed form of the primary cosmic ray 
flux in Eq. (|2|). Proton-proton collisions at these energies were then simulated using Pythia || 
resulting in complete events. After all particle decays, neutrinos and photons were recorded and 
filled in energy-histograms applying the weight for the event. 
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Figure 1: The E 3 -weighted flux of muon and electron neutrinos from interactions of cosmic rays 
with the interstellar medium in a column of length 1 kpc and density 1 nucleonj cm 3 . The total 
flux is shown (full curve) as well as the contributions from the decay of the indicated particles. 



3.1 The neutrino fluxes 

The resulting fluxes of muon and electron neutrinos are shown in Fig.|l|. For muon neutrinos the 
dominating source is muon and pion decay, while for electron neutrinos muon decay dominate 
strongly with only a few per cent from kaon decay and a totally negligible fraction from neutron 
decay. 

This is rather different from the other sources of cosmic ray induced neutrinos in the vicinity 
of the Earth, i.e. the Sun's and the Earth's atmosphere. In case of the Earth's atmosphere the 
distance involved is so short that the muon can be considered almost stable. The most important 
sources of muon (electron) neutrinos are therefore the decays of pions and kaons (kaons), with 
charmed particles being the dominating source at very high energies [Q]. A comparison with 
these atmospheric fluxes is made below in section 5 (Fig. 6). For the Sun, muons cannot be 
considered stable, and give a significant contribution to the fluxes. The muons, however, loses 
energy through electromagnetic interactions and are therefore less important at higher energies 
compared with the light mesons [14]. 

Because of the small cross section for z^-hadron and v— 7 interactions, the ^-fluxes will not be 
noticeably attenuated due to interactions with the interstellar medium or the cosmic microwave 
background. The attenuation due to interactions with stellar material is also negligible since, 
as mentioned, the stars cover only a very small solid angle. 

3.2 The photon flux 

As opposed to the neutrino fluxes, the photon or gamma flux is attenuated by interactions. Pho- 
tons of very high energy can interact with the cosmic microwave background to produce particles. 
Unless the energy is extremely high, one need only consider electron-positron production, i.e. 
77 — ► e + e~. 
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Figure 2: The attenuation factor, Eq. for the photon flux as a function of the photon energy 
Ky and the distance R from the Earth to the edge of the galaxy. 



The cosmic microwave background is isotropic and has an energy spectrum given by 

ME) = 1 £ 2 , Rno -5 _ F) 3-3 • 10 7 U) 

dE TT 2 (hc) 3 exp(E/kT) - 1 ^ { ' E ' 1 ' 

with E in eV . The first part is the thermal background spectrum at temperature T and the 
second part is a very simple parameterisation of the radio-wave background at very low energy 
[16]. The latter is far from as accurately know as the former, but since it only influences the 
flux above ~ 10 7 GeV it is not of significant importance for our purposes. 

To calculate the attenuation one needs to account for both the energy-spectrum of the 
background radiation and the angle between the two photons. We do this by defining an effective 
thickness (mbarn/cm 3 ) given by 

< an > (£L) = — P" deft P d{cos6) P dE^^-a^ fe (2E 1 E(l + cos6)), (5) 
Air Jo J-i JE t h(0) dE 

where (ft and 9 are the azimuthal and polar angles between the interacting gammas in the galactic 
'laboratory system'. is the cosmic background radiation in Eq. (|j) and c 77 ^^(s) is the 

cross section for charged lepton pair production as a function of the CM-energy y/s, which we 
have calculated. It is sufficient in this context to use the leading order formula for this cross 
section and to consider e + e~production only. The threshold energy for a microwave background 
photon is E t h(9) = 2mj/E~ / (l + cos(8)) in terms of the energy Ej of the incoming photon. 

Since the production of photons is isotropic along the line of sight within the plane of the 
Milky Way an attenuation factor can be derived and applied to the flux obtained from the Monte 
Carlo simulation. The attenuation factor is obtained by integrating over the line of sight and 
weighting each production point with the probability (e~ r<an> ) that a produced photon reaches 
the Earth, 
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Figure 3: The attenuated (solid line) and unattenuated (dashed line) E s -weighted flux of photons 
from interactions of cosmic rays with the interstellar medium of density 1 nucleon/ cm 3 integrated 
over the whole Milky Way in the direction to its center. 



where R is distance to the edge of the Milky Way. A numerical evaluation of this factor is 
shown in Fig. ||. The dip structure in the attenuation factor is a consequence of folding two 
peaked functions. The cosmic microwave background radiation which peaks at T ~ 7K, and 
the electron-positron production cross section peaked at s ~ Yim\. 

With this attenuation factor, we then obtain the attenuated photon flux shown in Fig. |3] and 
compared to the unattenuated flux; both are in the direction towards the centre of the Milky 
Way. The attenuation dip structure is obviously a reflection of the structure in Fig. 2. Another 
such dip structure would also appear at an energy of (m M /m e ) 2 10 6 ' 5 GeV ~ 10 11 GeV, when 
the reaction 77 — > becomes energetically possible. For such high energies is, however, our 

parameterisation of the primary flux inadequate and our method therefore no longer applicable. 

3.3 Parameterisation of the fluxes 

The resulting fluxes can be parameterised in a form inspired by the analytic formalism used 
by Berezinsky et al. [||. The photon flux, which is attenuated through pair production as 
discussed, is better treated in two steps; one for the production and one for the attenuation. 
The non-attenuated fluxes are parameterised as 

f RSNqE^- 1 E < E 
HE) = (7) 
[ RSN^E^'- 1 E > E 

with the fitted parameter values given in Table |l|. (Nq is not fitted, but given by the conti- 
nuity condition at Eq.) R5 is the column number density, 5 is the average number density in 
nucleons / err? and R is the distance to the edge of the galaxy in kpc. By inserting the corre- 
sponding values from our model, or another model, of the galaxy one can obtain the fluxes from 
this parameterisation. Results are shown in section || and discussed in connections with other 
sources. 
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Table 1: Values of the parameters in Eq. (fjp obtained from fits to the neutrino fluxes in Fig. [J 
and the unattenuated photon flux in Fig.\&. 



The attenuation factor A(E, R) can be obtained from Eq. (||) using the effective thickness 
parameterised as 

log{< na > (E )} = -25.34 + 6.82 Vlog E - 4.99 e~ a29 ( log s o-4.99) Eq > w 5 GeV ^ 

whereas for Eq < 10 5 GeV the pair production cross section negligible small and the attenuation 
factor effectively unity, see Fig. 0. 



4 Fluxes from a possible dark matter galactic halo 

Cosmic rays could also interact with matter in a galactic halo and, depending on the properties 
of this matter, give a neutrino flux of interest. If the halo matter consists of weakly interacting 
massive particles (WIMPS) the interaction cross section is very small and would give a negligible 
neutrino flux. However, with dark matter of hadronic nature this source could give a measurable 
flux. In the scenario of massive cold halo objects (MACHO's), the interactions would take place 
in objects of high density. The interaction lengths of the produced secondary mesons would then 
be shorter than their decay lengths such that they would loose energy or be absorbed before 
decaying into leptons. This would result in a lepton flux at very low energies only, which at 
the Earth would have to compete with the high atmospheric flux. Recent measurements by 
the MACHO collaboration [17] show, however, that these objects can only make up ~ 20% 



of the dark matter. This leaves the possibility for a gaseous hadronic component of the dark 
matter as discussed by De Paolis et al. H . They calculate gamma ray fluxes in the energy range 
1 — 10 6 GeV from decays of 7r°'s produced in cosmic ray interaction. 

Following the approach of De Paolis et al. || we have calculated the corresponding neutrino 
fluxes. The gaseous hadronic matter distribution in the halo is assumed to be of the form 

o 2 + R GC 



p(R) = Pof 



(9) 



+ R 2 ' 

where po ~ 0-3 nucleons / cm 3 is the local dark matter density and / ~ 0.5 is the fraction of dark 
matter in form of gaseous nucleons. R is the distance from the galactic centre to an arbitrary 
point and Rqc ~ 8.5 kpc to the Earth, whereas a ~ 0.5 kpc is the galactic core radius. 

The flux of cosmic rays in the halo may consist of two parts. First, particles that escapes 
from the disc, which probably give a flux that decreases in the same way as the halo density 
does, i.e. Eq. @, unless it is confined by some unknown mechanism. Secondly, there may be 
a flux of extra-galactic origin, which can be assumed to be isotropic. De Paolis et al. || do 
not consider energies higher than 10 6 GeV and are therefore dominated by the first part. The 
normalisation can then be estimated from theoretical arguments jyj resulting in ~ 1/500 times 
that at the Earth. Thus we have the cosmic ray flux in the halo 

,(£) a 2 + R 2 GC 



b p (E,R) 



500 



a 2 + R 2 ' 



(10) 
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Figure 4: The E 3 -weighted flux of muon neutrinos from a possible dark matter galactic halo 
induced by interactions of cosmic rays of galactic ( dashed line ) and extra- galactic ( dotted line ) 
origin. Shown for comparison is also the flux from the normal interstellar medium (section 3.1) 
at high galactic latitudes (solid line). 



Part of this flux is directed away from the galaxy such that the secondary flux will not reach 
the Earth. This could be accounted for by an efficiency factor, but since it is unknown it is not 
taken into account. 

The fluxes of neutrinos are then obtained by folding the cosmic ray flux and the halo density 
and integrating the production along the line of sight analogously to the galactic flux (Eq. ([!])). 
The resulting flux is shown by the dashed curve in Fig. [| and seen to be significantly smaller 
than the galactic flux at high latitudes (where it is lowest). One should remember that this 
flux from the halo is a crude estimate and may be seen as an upper limit due to the mentioned 
neglect of the efficiency factor. 

Extra-galactic cosmic rays are also believed to contribute to the flux in the halo, especially 
at very high energies. The normalisation of this component is unknown, but could be as high 
as given by the high energy flux at the Earth. In order to compare with the flux induced by 
cosmic rays escaping from the galaxy, we take the normalisation to be same as in Eq. (JToj) , 

M E,R)=^^1. (11) 

The neutrino flux is here obtained in the same way as above, and shown by the dotted curve in 
Fig.|. 

Although it seems that the flux from the halo cannot compete with that from the galactic 
disc, it is still conceivable that its magnitude is significantly higher than here estimated. This 
would happen if the primary cosmic rays are completely of extra-galactic origin and therefore 
the same as that at the Earth, i.e. Eq. ( pi] ) reduces to 4> P (E,R) = 4> P ,®{E). This would give a 
neutrino flux that is about a factor five higher than the corresponding interstellar flux. Although, 
such high extra-galactic flux is not probable, a higher flux than in Eq. ( pM| ) is not unrealistic. 

From the above discussion, the prospects for detection of a dark halo based on neutrino flux 
measurements does not look promising. In particular, this flux is considerably lower than the 
atmospheric flux [|J] and the estimated one from active galactic nuclei, as discussed in section 5 
(Fig. 6). 
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5 Concluding discussion 



The neutrino fluxes shown above in Fig.|l] are given for our basic 'unit column' of interstellar 
matter. Using these one may obtain the integrated flux from any direction. In Fig. [5] we show 
our resulting muon neutrino flux (full curves) from the direction towards the center of the Milky 
Way (highest flux) and orthogonal to the galactic plane (lowest flux). In Fig. ||a we compare with 
the original results of two earlier calculations by Domokos et al. |IJ and by Berezinsky et al. ^ . 
The differences between these curves are mainly related to the different assumptions concerning 
the interstellar matter density profile and the normalisation of the cosmic ray spectrum. Thus, 
they illustrate the uncertainty due to these inputs to the different calculations. 

In order show other differences between these different calculations we have recalculated the 
results of Q and [Q] using their formalisms but with our density and cosmic ray parametrisations. 
The results are shown Fig. ||b and demonstrate a close agreement for energies up to 10 5 -10® GeV . 
At higher energies there are, however, significant differences due to the treatment of the change in 
slope of the cosmic ray spectrum, cf. Eq. (|2|). In 0] energies above the 'knee' are not considered, 
and the naive extrapolation to higher energies reults in an overestimated flux. About a factor 
three excess at the highest energies is expected based on an estimate using the analytic method 
with spectrum- weighted Z-moments Q ||, [l9|], in agreement with the effect seen in Fig.||b. In |l|] 
this change of slope is included, but its effect on secondary particles cannot be fully taken into 
account since their calculation is based on an analytic method. One must here make assumptions 
about from what average primary energy a given neutrino comes, and thereby specify how the 
'knee' from the change in the cosmic ray energy spectrum is transported into a 'knee' in the 
final neutrino spectrum. This problem does not occur in our Monte Carlo method, since the 
neutrino spectra are here obtained through an event-by-event simulation taking the fluctuations 
into account. Our results are therefore more reliable in this respect. 
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Figure 5: The E 3 -weighted fluxes of interstellar muon neutrinos in the direction towards the 
center of the Milky Way (upper set of curves) and orthogonal to the galactic plane (lower set of 
curves). Comparison of our results (full lines) with those by Domokos et al.^j (dashed lines) 
and by Berezinsky et al.j$j (dotted lines). In (a) the original results are used, whereas in (b) the 
results of jfy and ||/ are modified to have the same galactic density profile and cosmic ray flux 
as in our calculation. 
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Figure 6: E 3 -weighted fluxes of muon neutrinos at the Earth from cosmic ray interactions with 
the interstellar medium looking towards the galactic centre (upper solid curve) and orthogonal to 
the galactic plane (lower solid curve ) as calculated in this study, and the flux from the Earth 's 
atmosphere in the vertical (V) and horizontal (H) ftsQ j direction originating from conventional 
it, K -meson decays as well as from prompt (P) charm decays For comparison, the diffuse 
flux of muon neutrinos from active galactic nuclei as predicted in fi2l\J (dotted line) and in J^/ 
(dash- dotted line) are shown. 



There is also a contribution to the interstellar flux from semileptonic decays of charmed and 
heavier hadrons. We have estimated this based on the analytic method with spectrum-weighted 
moments ||LJ| using our previoulsy calculated Z-moments for charm particle production 
and decay Taking also the contribution from muon decay into account, we find a muon 
neutrino flux which is only contributing about 2 • 10 -4 to the interstellar flux and a factor two 
higher contribution for the electron neutrino flux (due to a lower interstellar electron neutrino 
flux). The smallness of this charm contribution justifies the neglect of it in our Monte Carlo 
treatment. 

From Fig.||, and these considerations, one can also conclude that the uncertainties originat- 
ing from the particle physics input are not larger than other uncertainties. Our Monte Carlo 
calculation also confirmes that the approximation done with the analytic method are justfied 
except at the highest energies where special precautions must be taken. 

The fluxes from the interstellar medium should be compared with those from other sources 
in order to establish their observabelness and their significance as a primary interest of study 
or as a potential background. This is done in Fig. ^ where our result on the muon neutrino 
flux is compared with those from cosmic ray interactions in the Earth's atmosphere and with 
two predictions of the diffuse flux from active galactic nuclei. The interstellar flux is seen to be 
considerably lower than the atmospheric flux, except at the highest energies. It will therefore 
be very hard to observe with detectors on the Earth. At the highest energies considered, the 
interstellar flux is of comparable magnitude as the prompt atmospheric flux (from charm decays). 
The latter is almost direction independent up to ~ 10 7 GeV such that the curve shown for the 
vertical direction applies in essentially the whole energy range; the horizontal flux being slightly 
higher for the highest energies only. Given the small absolute scale of the fluxes at these energies, 
the event rate will be very low in the neutrino telescopes currently under construction. The 
interstellar flux becomes important compared to the vertical atmospheric flux only at energies 
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Figure 7: Energy-weighted spectrum of photons from interactions of cosmic ray particles, having 
the specified energies, with the interstellar matter (i.e. the photon yield). 



above 10 4 GeV . The corresponding event rate in a detector of 3 • 10 4 m 2 , e.g. Amanda, would 
be ~ 0.5 1 year in a cone of opening angle 10° directed towards the centre of the Milky way. 

Our calculated neutrino flux from the interstellar matter is significantly lower than the one 
from active galactic nuclei, as demonstrated in Fig.|(| and should therefore not be a problematic 
background for this extra-galactic source. On the other hand, the atmospheric neutrino flux will 
be problematic in this context, since it dominates in the lower energy range. At high energies 
where the atmospheric flux is lower, the absolute rate is very low. 

Concerning the flux of electron neutrinos the situation is very similar to the case of muon 
neutrinos, except that the atmospheric background is significantly lower Q and the interstel- 
lar flux is a factor ~ 2 lower than the corresponding muon neutrino flux (Fig. Since the 
dominating production process in AGN is decays of pions followed by the decays of the muons, 
the electron neutrino flux from AGN's should also be a factor ~ 2 lower. Thus, the relative 
fluxes from these sources should be the same as for muon neutrinos. There is, however, no good 
technique to detect high energy electron neutrinos, so the prospects to search for such a signal 
is worse. 

The yield of photons in cosmic ray interactions with the interstellar matter has been con- 
sidered before, e.g. by Berezinsky et al. [0]. Our corresponding result obtained from the unat- 
tenuated photon flux in section 3.2 is shown in Fig. 0. It is about a factor two higher than that 
of 0|, mainly due to the higher tt° multiplicity in our complete event simulations. The gamma 
flux can be used to set limits on the cosmic ray flux and the matter distribution in the galaxy. 
The experimental situation is here, however, quite different from the case with the neutrino flux. 
The gammas are detected either directly on satellite based experiments or with ground based air 
shower arrays. The signal rate will be higher than with neutrinos due to the larger interaction 
cross section. The flux of gammas is, however, very small (~ 10 -4 ) compared to the inclusive 
cosmic ray flux dominated by protons. Differences in the air showers produced by gammas and 
hadrons may allow access to information on the gamma flux. 
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